Quantum coherence and its dephasing in the giant spin Hall effect and nonlocal 
voltage generated by magnetotransport through multiterminal graphene bars 
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Motivated by the recent experimental observation [D. A. Abanin et a/., Science 323, 328 (2011)] of 
nonlocality in magnetotransport near the Dirac point in six-terminal graphene Hall bars, for a wide 
range of temperatures and magnetic fields, we develop a nonequilibrium Green function (NEGF) 
theory of this phenomenon. In the phase-coherent regime and strong magnetic field, we find large 
spin Hall (SH) conductance in four-terminal bridges, where the SH current is pure only at the Dirac 
point (DP), as well as the nonlocal voltage at a remote location in six-terminal bars where the direct 
and inverse SH effect operate at the same time. The "momentum-relaxing" dephasing reduces their 
values at the DP by two orders of magnitude while concurrently washing out any features away 
from the DP. Our theory is based on the Meir-Wingreen formula with dephasing introduced via 
phenomeno logical many-body self-energies, which is then linearized for multiterminal geometries to 
extract currents and voltages. 



PACS numbers: 72.80.Vp,72.15.Gd,72.10.Bg 

The recent experiments [1 on magnetotransport near 
the Dirac point (DP) in graphene have unveiled yet an- 
other exotic electronic property of this one-atom-thick 
carbon allotrope [2] which involves nonlocality and quan- 
tum mechanics while manifesting even at room temper- 
ature. The traditional observation of nonlocality, where 
current is injected through a pair of terminals and voltage 
is measured between another pair of terminals at some 
remote location, requires two-dimensional (2D) systems 
placed in a strong external magnetic field to generate the 
integer quantum Hall effect [3 (QHE) or spin-orbit cou- 
pling [4 (SOC) that can give rise to mesoscopic [5j |6] 
or quantum [3, 7 spin Hall effects (SHEs). In the for- 
mer case, nonlocality is due to transport through chiral 
edge states, while in the latter case injected longitudi- 
nal charge current generates transverse spin Hall current 
which is then detected in the remote part of the device 
via the inverse-SHE-induced voltages [6j[8] on the proviso 
that spins can survive dephasing between two locations. 

On the other hand, nonlocal voltage was observed in 
Ref. [1 even in weak magnetic fields B ~ 1 T and at 
room temperature T = 300 K which is outside of the 
integer QHE regime. Also, high mobility graphene sam- 
ples were supported by substrate made of atomically flat 
hexagonal boron-nitride that rules out Rashba SOC [4], 
introduced by charge impurities from the substrate [9] or 
lattice distortion by adatoms [10] , that would be respon- 
sible for the mesoscopic SHE scenario [5 . 

It turns out that SHE in the absence of SOC has a 
simple intuitive explanation based on classical Newto- 
nian dynamics of massless Dirac fermions. The classi- 
cal Hamiltonian of low-energy quasiparticles close to the 

Dirac point (DP) is given by H ± (p) = ±v F ^pl + P%, 
which in the weak external magnetic field B = V x A be- 
comes i^ ± (p) = ±v Fy J(p x - eA x ) 2 + (p y - eA y ) 2 . The 
classical velocity is then given by v^r = dH^/dpx = 




FIG. 1: (Color online) Schematic view of the six-terminal 
graphene Hall bar, modeled on the tight-binding lattice with 
single 7r orbital per site, which is employed to investigate non- 
local voltage between leads 5 and 6 due to current injected 
between leads 1 and 4. The dashed box on the left marks the 
four-terminal bridge used in the analysis of Zeeman-splitting- 
driven SHE where current injected between leads 1 and 4 
induces spin current in leads 2 and 3. The active region 
of the bar consists of graphene nanoribbon with armchair 
edges (AGNR) and a portion of semi-infinite leads modeled 
as GNRs with zigzag edges. For simplicity, external magnetic 
field, or many-body interactions responsible for dephasing, 
are present only within the illustrated active region. 

±v ^n^/Vn 2 , where II = p — eA, and the corresponding 
acceleration is 

, dv ev F v ± x B ev 2 7 v ± x B 

a = -r = / — ± ^ • (1) 

dt v/fP E± 

Thus, the quasiparticles with energy E + above the DP 
(or below with energy E~) moving in a weak (i.e., non- 
quantizing) perpendicular magnetic field will experience 
a transverse force which deflects them to the left (right). 
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FIG. 2: (Color online) The charge and spin transport quanti- 
ties in the four-terminal graphene bridge: (a) charge Hall con- 
ductance G xy = h/{Vi — V4); (b) charge Hall resistance Rn = 
(V 3 - V 2 )/h; (c) spin Hall conductance Gsu = 1$ /(Vi ~ V4); 
and (d) spin Hall angle 6> SH = ii/h- The width of AGNR 
channel is W/£b = 3.42 in the units of the magnetic length 
£b and a small "momentum-relaxing" dephasing d m = 0.047 
in introduced into the active region shown in Fig. [I] 



FIG. 3: (Color online) The charge and spin transport quan- 
tities in the four-terminal graphene bridge: (a) charge Hall 
conductance G xy ] (b) charge Hall resistance Rn] (c) spin Hall 
conductance Csh; and (d) spin Hall angle #sh- The width 
of the AGNR channel is W/£ B = 1.53 in the units of the 
magnetic length £b and large "momentum-relaxing" dephas- 
ing d m = O.47 in introduced into the active region shown in 
Fig.0 



Furthermore, when E ± is very close to the DP such de- 
flecting force will be very large. Although the Zeeman 
splitting Az in 2D electron gases (2DEGs) is typically 
small in a weak external magnetic field [4], it can play 
an essential role in graphene for /c#T < Az by shifting 
the Dirac cones for opposite spins to induce two types 
of carriers illustrated in the lower inset in Fig. [T] The 
quasiparticles with energy E + are spin-| polarized while 
those with energy E~ are spin-| polarized. These two 
effects, classical for charge and quantum for spin, con- 
spire to generate transverse spin current in response to 
longitudinal charge current, as illustrated in Fig. [I] Such 
phenomenology is similar to SHE in 2DEG bridges [5j [6] , 
even though no SOC is involved to provide the deflecting 
force of opposite direction for opposite spins [TT] . 

These simple arguments for the existence of the 
Zeeman-splitting-driven SHE (ZSHE) in graphene can 
be converted into a quasiclassical transport theory based 
on the Boltzmann equation [12]. However, quasiclassical 
theory is valid in high-T and weak- B regime, while ex- 
periments pQ have observed increasingly more profound 
nonlocality into the low-T and/or high-£? regime so that 
a unified theory is called for that can cover such wide 
range of parameters. For example, such theory should 
explain the nonlocal voltage in strong (quantizing) ex- 
ternal magnetic field but at intermediate temperatures 
where edge-state transport mechanism is removed. 

In this Letter, we develop a fully quantum transport 
theory of ZSHE in four-terminal graphene bridges, illus- 



trated by the device within a dashed box in Fig. [T] as well 
as the nonlocal voltage induced by the combination of di- 
rect and inverse ZSHE in six-terminal Hall bars shown in 
Fig. [l] This approach intrinsically accounts for the con- 
tributions of both electrons and holes, which is crucial 
to describe transport near the DP [12] . and it can also 
handle arbitrary scattering processes (in contrast to the 
Boltzmann equation which breaks down [13] close to the 
DP). Our central results, summarized in Figs.|2j[3]and[4j 
interpolate smoothly between the phase-coherent trans- 
port regime at low-T in the quantizing external mag- 
netic field and the semiclassical transport regime where 
dephasing by many-body interactions destroys features 
found at low-T while leaving peaks (of reduced magni- 
tude though) in the SH conductance and nonlocal volt- 
age around the DP in accord with experimental observa- 
tions pp. 

In the analysis of four-terminal bridges, voltage V/2 
is applied to lead 1 and — V/2 to lead 4 while voltages 
on leads 2 and 3 are set to zero. Figure [2] shows that in 
the quantizing external magnetic field W/£b > 1, where 
£b — y/h/\eB\ is the magnetic length in graphene, the 
four-terminal bridge generates large spin Hall conduc- 
tance Gsh = (^2 — ^sb/O^i — ^i) an d the corresponding 
SH angle # S h = G S h/G l with G L = V(Vi - V 4 ) being 
the longitudinal charge conductance. The spin current 
I2 = I2 — I2 in the ZSHE is carried by spins polarized 
along the z-axis orthogonal to the plane of graphene. The 
value of Gsh is comparable to the one predicted [14] for 
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FIG. 4: (Color online) Panels (a) and (c) plot charge Hall 
resistance Ru = (Vi — V4V/3, while panels (b) and (d) plot 
nonlocal resistance -Rnl = (V5 — Vq)/I\ as the central quan- 
tity measured in the recent experiments pQ on six-terminal 
graphene Hall bars. The quantum coherence is retained in 
panels (a) and (b) where only a small "momentum-relaxing" 
dephasing d m — O.O27 is present in the active region of the 
bar, while much larger dephasing d m = O.57 is used for pan- 
els (c) and (d). The width of the AGNR channel in Fig. [I] 
is W/£ B = 3.42 in panels (a) and (b) and W/£ B = 1.53 in 
panels (c) and (d) in the units of the magnetic length £b- 

the Rashba SO-coupled four-terminal 2DEGs of the size 
of the spin precession length (on which spin-^ precesses 
to spin-^ state). Unlike mesoscopic SHE [5 in 2DEGs 
where Rashba SOC induces both the transverse spin de- 
flection [11] and spin dephasing which compete against 
each other in the processes of generating pure (not accom- 
panied by any net charge flux) spin current, in the ZSHE 
transverse spin current is pure only at the DP [where 
charge current I2 = l\ + ^2 becomes zero in Figs. |i|a) 
and [3^a)] and spin precession is absent This might be 
advantageous for spintronic applications since spin de- 
phasing is evaded, as demonstrated by the experimental 
detection of nonlocal voltage even at distances ~ 10 jam 
away from the device region where SH current was in- 
duced p.!. We note that for very strong magnetic field, 
as could be achieved in ferromagnetic graphene, the Gsh 
peaks in Fig{2|c) would become quantized [15 as a real- 
ization of quantum SHE |3] in the absence of SOC. 

The introduction of dephasing processes into the four- 
terminal bridge, which relax both [16] the phase and the 
momentum of quasiparticles propagating through the ac- 
tive region, destroys the quantization of the charge Hall 
conductance G xy = h/(Vi — V4) or charge Hall resistance 
i?H and underlying chiral edge states, as demonstrated by 
the transition from Fig.[2|a) to Fig.|3|a) for G xy and from 
Fig. [2^b) to Fig. [3jb) for R^. The charge Hall resistance 
in four- terminal bridges is defined as Ru = (V3 — V2j/I\ 



for the measuring setup where current Ii is injected into 
lead 1 and voltages V3 and V2 develop as the response to 
it. The SH conductance and SH angle are concurrently 
reduced by two orders of magnitude, which are values 
similar to those found in quasiclassical approaches [12] 
in the temperature range T = 200-300 K. 

In the analysis of six-terminal Hall bars, charge cur- 
rent Ii is injected through lead 1 and current — 1\ flows 
through lead 4 while I a = in all other leads. We then 
compute voltages which develop in the leads a = 2, 3, 5, 6 
labeled in Fig. [I] in response to injected current 1\. Fig- 
ure [4jb) shows peaks in the nonlocal resistance, defined 
as .Rnl = (^5 — Vq)/Ii, within the phase-coherent trans- 
port regime which closely resemble the DP and side peaks 
observed experimentally in strong (quantizing) external 
magnetic field [1 . The transition of Rnl from Fig. |4^b) 
to Fig. |4jd) shows how dephasing removes both side 
peaks while leaving the nonlocal voltage around the DP 
which is two orders of magnitude smaller than in the 
phase-coherent regime. The Hall resistance in the six- 
terminal bar, i?H = (Vi — V^/Is defined for current in- 
jected Is and voltages measured between leads 1 and 4 
(for ii = I4 = 0), changes smoothly from Fig. [4ja) to 
Fig.Kc) as dephasing in increased, where the curve in 
Fig. |4|c) looks exactly the same as those observed exper- 
imentally for T = 250 K and B = 1-12 T [12]. 

In the rest of the paper, we discuss details of the mi- 
croscopic Hamiltonian model for multiterminal graphene 
Hall bar in Fig. [T] and the nonequilibrium Green function 
(NEGF) formalism including dephasing which is applied 
to understand transport properties of this device. Close 
to the DP, graphene can be can be described by the tight- 
binding Hamiltonian with single 7r-orbital per site 

n (nm), a 

(2) 

Here e n is the on-site energy, a = +1 for spin-^ elec- 
tron and <j = —1 for spin-| electron so that Zeeman 
splitting is given by Az = 2g/i j e£?, cj^ (c ncr ) creates 
(annihilates) electron with spin a in the 7r-orbital lo- 
cated on site n, and 7 is the nearest-neighbor hopping 
parameter. The external magnetic field enters through 
the phase factor (j) nm = ^ f ™ A • ds where the vector 
potential A = (By, 0, 0) is chosen in the Landau gauge 
and </>o = h/e is the flux quantum. The weak vs. strong 
magnetic field is tuned using the ratio W/£b, where W 
is the width of the AGNR channel of the bar in Fig. [T] 
All graphene bars studied in Figs.[2j [3] and [4] are placed 
in the quantizing magnetic field, W/£b > 1- 

While the integer QHE and quantum SHE have in- 
troduced the intricate physics of topologically ordered 
phases [3], their operational description used to analyze 
transport measurements [6l [7] is typically based on the 
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multiterminal Landauer-Biittiker (LB) formulas [17] 



i a = jJ2 T <*?( v <*- v ^> 



(3) 



written here assuming zero temperature. They relate 
charge current in lead a to voltages Vp in all other leads 
attached to the sample via the matrix of transmission 
coefficients T a p. These formulas are valid when phase 
coherence is retained in the active region of the device, 
while phase breaking events are assumed to be taking 
place only in the reservoirs to which the leads are at- 
tached at infinity and where electrons are equilibrated to 
acquire the Fermi-Dirac distribution f a (E) = f(E—eV a ). 

To take into account dephasing effects phenomenologi- 
cally, Biittiker introduced [18 an elegant concept of volt- 
age probes attached to the active region where no net 
current flows through them, so that for every electron 
that enters the probe and is absorbed by its reservoir 
another one has to come out which is not coherent with 
the one going in. For example, to apply this method 
to the graphene Hall bar in Fig. [I] one can attach one- 
dimensional leads to each site [19 of the honeycomb lat- 
tice. This is equivalent to adding a complex energy — ir] 
to e n in the Hamiltonian Eq. Q (parameter rj is related 
to the dephasing time r] = H/2r ( f } ). In addition, one has to 
solve Eq. (|3| by imposing that current through extra ID 
leads is zero thereby ensuring conservation of the total 
charge current. 

However, besides washing out quantum-coherence- 
generated fluctuations in T a p, Biittiker voltage probes 
are also introducing additional scattering (i.e., reduc- 
tion of T a p) in an uncontrolled fashion [16]. The NEGF 
formalism [20] provides a rigorous prescription for in- 
cluding any dephasing process to any order by starting 
from a microscopic Hamiltonian and by constructing in- 
teracting self-energies due to electron-electron, electron- 
phonon [20] or electron-spin [21] interactions. Although 
the NEGF formalism is virtually the only fully quan- 
titative quantum transport approach capable of scaling 
to large systems [22], the self-consistent computation of 
self-energies by starting from some microscopic many- 
body Hamiltonian is at present prohibitively expensive 
for devices containing realistic number of atoms. Thus, 
to include dephasing processes in the device in Fig. [I] 
containing few thousands of carbon atoms, we adopt a 
phenomenological model of Ref. [16] that is comparable 
to Biittiker voltage probes in conceptual and numerical 
simplicity, and yet allows one the flexibility of adjusting 
the degree of phase and momentum relaxation indepen- 
dently. 

The two fundamental objects [20] of 
the NEGF formalism are the retarded 

Qm'tM') = ~i®(t - *0({CnaW,4a'(*')}) ^d the 

lesser G<£ a ' (t, t') = ^^'(O^M) GFs which de- 
scribe the density of available quantum states and how 



electrons occupy those states, respectively. Here (...) 
denotes the nonequihbrium statistical average [20] . 
In stationary problems, G r and G < depend only on 
the time difference t — t' or energy E after Fourier 
transformation. Their matrix representations in the 
basis of local orbitals, such as the tt ones in Eq. (|2|, 
satisfy the following equations 



G r (E) 



£-H-££;(£-eK)-SL t (£) 



(4) 



G<(£) = G r (E) 



G a (E). (5) 



Here Yi r a (E) is the retarded self-energy determining 
the escape rates for electrons to exit into the attached 
leads, E<(25) = -f a \ET a {E - eV a ) - - eV a )\ 

is the corresponding lesser self-energy matrix due to 
the coupling to the leads, and advanced quantities 
are defined by O a = [O 7 "]**". In the "momentum- 
conserving" model of dephasing, the interacting 
self-energies are given by Yi\ nt (E) = d p G r (E) and 
T,^ t (E) = d p G < (E) ^ while in the "momentum- 
relaxing" model EFntCE) = WmG r (E)) and 
££ t CE) = V[d m G < (E)} [16. The operator £>[...] 
selects the diagonal elements of the matrix on which it 
acts while setting to zero all the off-diagonal elements. 
Any linear combination of these two choices can be 
used to adjust the phase and momentum relaxation 
lengths independently. When computed self-consistently 
together with G r (E) and G < (E), both of these choices 
for J}T nt (E) and J^^ t (E) ensure the conservation of 
charge current, ^2 a I a = 0. 

The "momentum-relaxing" model of dephasing we 
adopt here accounts for local simultaneous phase and mo- 
mentum relaxation, and it can be interpreted as a highly 
simplified version (valid in the high-temperature limit) of 
the so-called self-consistent Born approximation [20] for 
electron-phonon interaction. This model has also been 
employed before to study dephasing effects in the integer 
QHE [23] where phenomenological dephasing length is 
often employed [24] to account for electron-electron and 
electron-phonon interactions without delving into micro- 
scopic details. 

Assuming that dephasing is localized within the ac- 
tive region of the graphene Hall bar, the spin-resolved 
charge current in lead a is given by the Meir-Wingreen 
formula [20] 

1% = | JdETr[^ aa (E)G > ^ a (E)-^ aa (E)G < ^ (T (E)}. 

(6) 

The total charge current in lead a is I a = 1^ + 1^ and 
the total spin current is if = l£ — /J;. The first term 
in Eq. (|6| gives the current flowing from lead a towards 
the active region (because it is proportional to G y (E) 
which describes the empty states in the active region), 
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while the second term gives the current flowing in the 
opposite direction (because it is proportional to G < (E) 
which describes the occupied states in the active region) . 
Likewise, the self-energies E^(i2) are proportional to the 
occupied lead states and the empty lead states, respec- 
tively. 

While Eq. |6| is valid both in the linear and non-linear 
transport regimes, in its original form it is not useful 
for the analysis of currents and voltages in multitermi- 
nal devices. That is, instead of voltages hidden in the 
self-energy and GF matrices, one would like to recast 
Eq. (|6| into the form similar to Eq. ([3| where one can 
easily invert such equations to obtain voltages measured 
between the terminals for known currents injected into 
the device. For this purpose, we expand all quantities 
in Eq. Q and (|5| to linear order in V a \ {%) G r (E) w 

G5(£0 + G5(£7)E a - eV a ) - V r a (E)]G r (E); (»0 

17 a {E - eV a ) « 17 a {E) - eV a d!7 a {E) / dE- and (iii) 
f a (E) « f(E) - eV a df/dE. Here G r (E) = [E - H - 
X)^ (£?)]" 1 is the retarded GF in equilibrium, V a = const. 
Using this in Eq. (|6|, together with the expressions for 
J2l nt (E) and 5]^ t (£') for "momentum-relaxing" dephas- 
ing, yields the following generalization of Eq. ([3| 

7 « = T E( T « / + T aP° h )(Va ~ Vp). (7) 

The coherent transmission coefficient is = 
Tr {r^Gor^Go}, and the incoherent contribution is 
T^ coh = Tr{r a G5r^Gg}. Using the notation [A] jv 
for the matrix element of A, the diagonal elements 
[Tp]jj = d m J2 v [Q] jv [G^TpG^] vv are expressed [23] in 
terms of Q = [1 - d^P]' 1 and [P] jv = [G r ] jv [G%] vj . 

In conclusion, we have developed a fully quantum 
transport theory of recently observed [1 nonlocal volt- 
age in magnetotransport near the DP in graphene Hall 
bars which provides a unified picture of this phenomenon 
and the underlying ZSHE from the phase-coherent trans- 
port regime at low temperatures to semiclassical trans- 
port regime at higher temperatures while allowing one 
to take into account arbitrary strength of magnetic field 
or scattering processes near the DP. Our theory starts 
from the NEGF-based Meir-Wingreen formula, including 
phenomenological many-body self-energies that take into 
account relaxation of both the phase and the momen- 
tum of Dirac fermions in the active region of the device, 
which is then linearized to provide connection between 
current and voltages in different leads thereby generaliz- 
ing the usual LB formulas for phase-coherent transport 



in multiterminal geometries. 
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